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ESTIMATION  OF  TOE  MECHANICAL  PROPERTIES  OF  ROCK  MASSES 


UTOODUCTION 

The  Review  of  Progress  of  Contract  No.  H0101610,  ’’The  Mechanical 
Properties  of  Rock  Masses",  submitted  on  May  28,  1971,  [two  weeks  prior 
to  the  end  of  the  Contract  period  (June  11,  1971)3  intimated  (Note 
on  Page  i)  that  the  Annual  Technical  Report  would  contain  "an  extensive 
literature  review  and  more  detailed  information  on  all  phases  of  the 
research" . 

After  discussion  with  the  A.R.P.A.  Project  Engineer  for  this  contract, 
it  was  decided  that  the  Final  Report*  should  instead  consist  of  a  concise 
general  review  of  the  overall  aims  and  achievement?  of  the  one  year  program 
of  research.  This  is  the  purpose  of  the  remarks  below. 

i'liC'l  IANICAL  PROPERTIES  OF  ROCK  MASSES 

Rational  design  of  engineering  strictures  in  or  on  rock,  and  the  under¬ 
standing  of  rock  mass  response  to  dynamically  and  statically  applied  loads s 
requires  a  knowledge  of  the  mechanical  properties  of  the  rock  mass  involved. 
Reliable  specification  of  these  properties  is  very  difficult  to  achieve, 
for  several  reasons : 

'hie  size  of  the  mass  for  which  the  properties  need  to  be  known  is  frequently 
>'f  the  order  of  hundreds  of  feet,  or  larger.  To  test  such  a  region  over  the 
range  of  loads  to  which  the  prototype  will  be  subjected  is  not  feasible.  Apart 
from  the  astronomical  cost  that  would  be  involved,  the  test  itself  could  cause 

i - 

This  replaces  the  Annual  Technical  Report,  since  no  funds  were  awarded 
for  continuation  of  Contract  No.  H0101610  for  the  second  year  of  the 
proposed  work. 
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peimanent  damage  of  the  mass .  Since  no  two  geological  sites  are  identical 
it  would  be  Inadviseable  to  test  any  mass  other  than  that  actually  involved 
In  the  engineering  problem. 

The  only  realistic  approach  is  to  make  sarrple  determinations  of  the 
needed  properties,  and  then  make  predictions  of  the  mass  response  to  the 
loads  applied  In  practice.  The  actual  behaviour  should  then  be  monitored 
by  instruments  installed  according  to  a  plan,  derived  together  with  the 
predictions,  to  check  the  response  of  the  mass  against  that  predicted. 

In  this  way,  the  analytical  procedures  can  be  refined,  and  progress  made 
towards  the  reliable  design  of  engineering  structures  in  rock.  It  is 
generally  accepted  that  che  mechanical  behaviour  of  a  rock  mass  is  deter- 
mined  by  the  properties  of  its  two  main  structural  components,  viz. 

1.  relatively  intact  [volumetric]  rock  elements; 

2.  structural  discontinuities  [more  or  less  planar]  which  bo'^nd 
the  elements. 

The  mechanical  properties  involved  In  practical  problems  also  fall  into 
two  broad  categories,  concerned  with  either 

1.  deformability ,  or 

2.  strength 

of  the  structural  components. 

Several  methods  have  been  developed,  both  in  the  field  and  in  the 
laboratory,  intended  to  measure  these  properties. 

Field  tests  are  valuable  in  that  they  test  the  rock  'in-situ'  under  more 
or  less  the  same  environmental  conditions  as  the  mass  [although  significant 
changes  (e.g.  In  stresses  acting  on  the  test -piece)  may  occur  during  develop- 
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ment  of  the  site!].  It  is  also  possible  to  conduct  tests  on  a  somewhat 
larger  scale  than  is  usual  in  the  laboratory.  Such  field  tests  are 
correspondingly  expensive.  One  radial  jacking  test,  for  example,  in 
which  an  8  ft.  length  of  an  8  ft.  diameter  tunnel  is  loaded,  may  cost 
around  $50,000  -  $75,000.  Determination  of  the  strength  of  an  under¬ 
ground  coal  pillar  requires  comparable  expenditures. 

Laboratory  measurements  allow  tests  under  more  closely  controlled 
loading  conditions,  and  can  be  repeated  to  permit  estimation  of  the  re¬ 
producibility  of  data.  The  mechanical  behaviour  of  the  specimen  can  be 
explored  over  a  wide  range  of  loading  conditons;  and  costs  are  relatively 
low.  Current  understanding  of  the  mechanics  of  rock  fracture  is  quite 
poor.  For  example  Ring  tests  and  Brazilian  tests  are  ccnmonly  used  to 
determine  the  tensile  strength  of  rock  specimens.  Values  computed  in  the 
standard  manner  for  each  method  convnonly  differ  by  300%  or  more  for  the 
same  rock  type.  'Hie  same  is  true  for  compression  strength.  This  fact  is 
Less  well  recognized  since  congress  ion  test  procedures  are  more  standardized 
than  tension  tests. 

In  all  cases  (i.e.  both  field  and  laboratory)  measurements  are  made  on 
samples  which  are  usually  quite  small  compared  to  the  prototype,  and  serious 
questions  are  raised  as  to 

(a)  the  correlation  between  the  sample  results  and  the  corresponding 
properties  of  the  full-scale  conponent,  and 

(b)  the  relationship  between  the  full-scale  conponent  properties  and 
the  overall  properties  of  the  rock  masses. 


Engineers  assume  that  strength  decreases  and  defonnability  increases 
with  increase  in  size  of  the  loaded  rock  mass,  but  there  is  little  under¬ 
standing  of  the  extent  of  the  ’size*  effect.  There  is  conclusive  evidence 
that  the  strength  of  coal  pillars  drops  appreciably  with  size  increase, 
but  coal  is  not  typical  of  sedimentary  or  igneous  rocks. 

Unless  such  questions  are  answered,  then  rock  testing  results  will  continue 
to  be  treated  as  little  more  than  rough  qualitative  indices.  Current  move:-, 
to  standardise  laboratory  testing  procedures  (welcome  in  one  sense)  tend  to 
promote  this  attitude.  Much  laboratory  research  in  rock  mechanics  has  de¬ 
generated  into  a  striving  for  2?  reproducibility  of  results  for  a  particular 
test  procedure,  with  apparent  disregard  for  the  fact  that  discrepancies  of  the 
order  of  300?  remain  unexplained.  Such  differences  indicate  that  the  mechanics 
of  rock  fracture  are  not  understood.  Until  they  are  the  significance  of  the 
numerical  value,  on  which  so  much  effort  is  expended,  is  unknown.  Other  load¬ 
ing  conditions  cannot  be  interpreted  in  terms  of  such  tests.* 

In  the  field  of  soil  mechanics  which  is  much  more  older  than  rock  mechanics , 
Professor  Roscoe  has  stated  (Tenth  Rankine  Lecture,  Geotechnique ,  20,  No.  2, 

Page  129),  with  reference  to  the  aims  of  the  Cambridge  soil  mechanics  group, 
that,  "The  principal  objective  is  to  develop  an  understanding  of  the  stress -strut 
behavior  of  soils  so  that  reliable  predictions  can  be  made  concerning  their 
load-deformation  characteristics  at  all  working  loads..."  and  it  is  considered 
that  this  philosophy  is  equally  applicable  to  rock  mechanics;  the  deformation  and 
failure  behavior  of  rock  must  be  studied  in  totality  rather  than  with  respect 
to  specific  characteristics  of  isolated  tests. 

•Attempts  liave  been  made,  for  exanple,  to  relate  the  energy  required  fx*r  unit 
volume  of  rock  broken  in  drilling  (E^)  to  the  rock  'compressive  strength’  (uc). 
Results  indicate  that  0. 3oc<Ev<3.0oc,  i.e.  that  there  is  1000?  variation  In  the 
ratio  By/oc!  This  exanple  also  illustrates  the  extent  to  which  the  compressive 
strength  is  looked  upon  as  a  basic  property .  Ey  can  be  obtained  directly  in 
the  laboratory,  simply  by  drilling  a  sample  of  the  rock.  Such  a  test  is  quicker, 
easier,  and  obviously  a  much  more  reliable  index  of  drilling  performance.  The 
compressive  strength  value  of  the  rock  is  invariably  requested. 
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It  was  with  these  questions  In  mind  that  a  proposal  was  submitted  to  A.R.P..A. , 
suggesting  that  research  was  needed  to  determine 

(a)  the  intrinsic*  mechanical  properties  of  a  rock  (i.e.  properties 
whose  magnitude  is  unaffected  by  experimental  technique) 

(b)  rational  procedures  for  estimating  the  mechanical  properties 
of  rock  masses  from  known  properties  of  the  structural 
components  of  the  mass. 

Although  Contract  H0101610,  which  resulted  from  the  proposal,  requested 
some  particular  experiments,  it  did  in  general,  request  a  research  effort 
aimed  towards  the  determination  of  (a)  and  (b)  above  over  the  two-year  period 
June  12,  1970  -  July  11,  1972. 

The  research  surrmarized  below  indicates  what  was  accomplished  towards  these 
goals  in  the  first  year  of  the  proposed  program.  It  is  important  to  note  that 
the  research  was  planned  with  a  two  year  program  in  mind,  so  that  some  of  the 
work  liad  not  reached  any  definite  conclusions  by  the  end  of  the  first  year. 

INTRINSIC  PROPERTIES  OF  ROCKS 

We  will  consider  the  mechanical  properties  in  the  two  categories  already 
mentioned,  viz.  deformation  properties,  and  strength  properties.  It  is  con¬ 
venient  to  consider  strength  properties  first. 

STRENGTH  OF  ROCK 

The  most  rational  and  physically  attractive  explanation  of  the  strength 
of  brittle  materials  [a  class  into  which  most  rocks  can  be  placed,  at  least 

as  an  approximation]  is  that  first  offerred  by  Griffith  in  1921.  Although  the 

* - 

As  distinct  from  'experimental  properties'  such  as  tensile  and  conpresslve 
strength,  values  which  appear  to  vary  considerably  depending  on  the  test 
procedure  used. 
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Griffith  theory  has  since  been  found  to  be  an  over-simplification  of  the 
pi  jy  si  cal  processes  of  fracture  the  central  notion  on  which  it  is  based, 
i.e.  that  fracture  in  a  material  develops  with  the  onset  of  unstable 
energy  changes  in  the  material,  is  of  very  general  validity.  It  is 
essentially  no  more  or  less  than  the  application  of  the  Theorem  of  Minimum 
Potential  Energy  to  fracture. 

This  theorem  applies  to  rock  fracture  as  much  as  to  other  physical 

processes.  It  is  important  to  note  that  it  applies  on  any  scale.  If 

it  can  be  used  to  describe  fracture  on  the  small  scales  it  can  be  used  to 

describe  fracture  on  the  large  scale  Griffith  was  particularly  concerned 

with  the  problem  of  strength  on  the  small  scale  and  so  restricted  his 

attention  to  that  scale.  In  applying  the  Theorem  of  Minimum  Energy  to 

fracture  he  states  the  following  principle: 

...the  equilibrium  position,  if  equilibrium 
is  possible,  must  be  one  in  which  rupture 
of  the  solid  has  occurred,  if  the  system 
can  pass  from  the  unbroken  to  the  broken 
condition  by  a  process  involving  a  continuous 
decrease  in  potential  energy. 

There  is  no  mention  of  size  in  '^his ,  very  general ,  statement . 

The  Principal  Investigator  has  in  fact  outlined  how  the  principle  may 
be  applied  to  explain  the  reduction  in  ’strength'  with  size  of  underground 
tunnels  and  how  it  may  be  employed  in  the  prediction  of  fracturing  an!  cratering 
in  blasting. 

Of  specific  relevance  to  Contract  No.  H0101610  the  principle  implies  also 
that  if  the  energy  instability  approach  can  be  shown  to  be  of  practical 
value  in  explaining  laboratory  strength  results,  it  should  be  equally  applicable 
in  the  analysis  of  rock  mass  stability  and  strength. 


-7- 


Application  of  energy  Instability  analysis  to  fracture  requires  a 
knowledge  of  the  characteristic,  or  intrinsic,  mechanical  properties  of 
the  material  which  control  the  energy  demand  -  energy  supply  relationship. 

The  search  for  intrinsic  material  properties  of  rock  is  thus  synonymous 
with  a  test  of  the  practical  utility  of  the  energy  stability  analysis  in 
explaining  rock  fracture. 

It  is  in  the  context  of  the  statements  in  the  last  two  paragraphs  that 
the  significance  of  the  laboratory  tests  in  Contract  H0101610  should  be  judged. 

STRENGTH  OF  ROCK  SPECIMENS 

To  begin  the  test  of  the  theory  on  the  laboratory  scale  is  appropriate 
for  several  reasons.  Thus, 

1.  The  theory  is  general  and  applicable  on  all  scales.  Experimental 
proof  of  validity  on  one  scale  is  evidence  of  general  validity. 

2.  Laboratory  tests  are  relatively  inexpensive,  and  the  theory  can 
be  checked  for  a  wide  variety  of  loading  conditions. 

3.  There  is  a  wealth  of  laboratory  data  already  available  which,  when 
analyzed  by  existing  [stress  limit]  criteria  of  rock  failure,  are 
quite  inconsistent.  Hie  value  of  the  energy  instability  theory  can 
be  quickly  assessed  by  using  it  to  re-interpret  these  data. 

'Hie  specific  laboratory  studies  were  selected  to  check  the  predictions  of  the 
Griffith  theory,  expanded  and  reinterpreted  so  as  to  be  more  directly  applicable 
to  laboratory  measured  values. 

The  literature  abounds  in  detailed  reviews  of  the  Griffith  theory  so  that 
it  would  be  superfluous  to  present  an  extensive  treatment  of  it  here.  It  is, 
however,  instructive  to  review  the  main  results  and  implications  for  laboratory 


testing  and  rock  fracture. 


Griffith  Theory  of  Tensile  Rupture 

The  observed  tensile  strength  of  brittle  materials  is  generally  one, 

two  (or  even  more)  orders  of  magnitude  lower  than  the  theoretical  inter- 

molecular  cohesive  strength.  The  latter  may  be  shown  to  be  of  the  order 

of  E/10,  where  E  is  the  modulus  of  elasticity  in  tension  for  the  material. 

Most  granites,  for  example,  have  a  tensile 
modulus  around  7  x  10°  p.s.i.  This  would 
Indicate  a  theoretical  tensile  strength  of 
around  700,000  p.s.i.  Measured  'tensile 
strength'  is  of  the  order  of  2,000  p.s.i., 
or  less  than  1/300  of  the  theoretical 
value . 

Griffith  explained  this  discrepancy  as  due  to  the  existence  of  flaws  in 
the  material.  These  flaws  served  to  concentrate  the  applied  load  to  very 
high  values  (stresses)  at  tlie  crack  tips.  The  high  stresses  were  necessary 
to  overcome  the  intermolecular  attraction  [cohesive  strength]  but  alone  were 
not  sufficient  to  cause  crack  extension  (i.e,  fracture).  For  this  to  occur, 
the  potential  energy  of  the  system  [i.e.  potential  energy  of  the  applied  load: 
plus  stored  (strain)  energy  of  the  loaded  body,  plus  energy  needed  to  create 
the  new  surface  of  the  extended  crack]  must  decrease  continuously  with  crack 
extension. 

The  second  requirement,  i.e.  the  energy  condition  is  the  more  important 
for  rock,  since  the  necessary  stress  concentrations  are  undoubtedly  present 
due  to  the  (small-scale)  inhomogeneous  structure  of  all  rocks. 

The  (plane  stress)  situation  analyzed  by  Griffith  is  shown  in  Figure  1. 

A  sharp  crack,  length  2c,  radius  of  curvature  p  at  the  tip,  is  introduced 


o 


into  an  infinitely  large  plate,  [later  we  shall  assume  the  plate  to  have  a 

large  but  finite  planar  area  A]  of  a  homogeneous  isotropically  elastic., 

material.  The  plate  has  unit  tliickness  and  is  uniformly  stressed  at  the 

outer  boundaries  AB,  CD,  to  an  intensity  'o',  under  dead-load  conditions 

(i.e.  the  stress  does  not  change  with  displacement  of  the  loaded  boundaries  - 

tills  will  later  be  referred  to  as  an  'infinitely  soft'  loading  system). 

'Die  change  In  potential  energy  of  the  system  (AP)  which  occurs  due 

to  introduction  of  the  crack  is  found  to  be 

AP  =  ^Icy  -  ttc  o  (1) 

E 

Tiie  first  term,  4cy,  is  the  increase  in  surface  energy  due  to  introduction  of 

tiie  2  crack  surfaces  each  2c  in  area  (y  is  the  specific  surface  energy).  The 
2  2 

second  term,  -ttc  o  ,  is  the  reduction  in  the  potential  energy  of  the  remainder 
E 

of  the  system  [i.e.  change  in  the  potential  energy  of  the  applied  forces,  plus 
change  in  strain  energy  of  the  loaded  plate]  due  to  introduction  of  the  crack. 
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2  2 

It  may  be  shown  [see  publication  (4)#]  that  this  quantity,  -vc  o"  ,  is 

E 

constant,  independent  of  the  stiffness  of  the  applied  load.  [This 

'stiflYiess'  is  the  force-displacement  relationship  for  the  applied 

load.]  The  relationship  of  Eq.  1  is  shown  graphically  for  two  levels 

(o  ,  o,  )  of  applied  stress,  with  a  <  a.  [tension  assumed  positive), 
an  a  d 


Figure  2 .  Change  in  potential  energy  of  a  cracked  plate 
[as  in  Figure  1]  with  crack  extension;  for  two 
levels  a  ,  o^,  of  applied  stress. 


S - 

References  so  noted  indicate  publications  resulting  wholly  or  partially 
from  Contract  H0101610  support.  Hie  list  of  contract  publications  is 
given  in  Appendix  1. 
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According  to  the  Griffith  criterion,  fracture  will  occur  [i.e.  the  crack 

will  extend]  only  when  the  potential  energy  of  the  system  decreases  with 

increase  in  crack  length,  i.e. 

3_(AP)  <  0  (2) 

3c 

Thus ,  referring  to  Figure  2,  for  constant  applied  stress,  a  the  crack  will 

a 

not  extend  for  crack  lengths  less  than  2c  .  Similarly  at  the  higher  stress 

a(j  the  crack  will  not  extend  for  cracks  less  than  2c^  in  length. 

Applying  the  criterion.  Equation  (2),  to  the  energy  change  expression 

squat  ion  (1),  we  find 

3  (AP)  =  -  2ttcc2  <  0 

3c  E 

(3) 

or  a  -  1 2Ey 

Jnc 

o  is  the  stress  level  at  which  spontaneous  extension  of  a  crack,  original 
length  2c,  will  occur,  i.e.  the  plate  will  fracture. 

The  above  analysis,  has  been  made  for  plane  stress,  but  vet 7  similar 
results  are  obtained  for  other  loading  condition;,  and  flaw  shapes.  The 
important  point  is  that  all  are  of  the  fonn 

a^c  =  const.  (;l) 

which  is  inherent  in  the  type  of  elastic  Instability  analysis. 

The  Griffith  theory  has  been  objected  to  on  the  basis  that  the  predicted 
value  of  the  tensile  strength  obtained  by  insert ing  'appropriate'  values  ot 
surface  energy'  (y),  modulus  of  elasticity  (E)  and  semi -crack-length  (c) 
into  liquation  (3),  are  too  low.  It  is  now  well  recognized,  for  metals,  at  least , 
that  the  energy  dissipated  in  crack  extension  is.  much  greater  than  'y'  ;  a 
considerable  amount  of  plastic  work  is  involved. 
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Tlie  important  point  to  draw  out  of  the  Griffith  theory  is  the  following 

fracture  will  occur  when  the  rate  at  which  energy  is  released  by  the 
loaded  system  exceeds  the  rate  at  which  it  is  dissipated  by  the  processes 
of  fracture  extension. 

We  may  generalize  Equation  (1)  to  the  form 

AP  =  W '  -  W  ('•) 

d  s 

where  AP  is,  as  before,  the  increase  in  potential  energy  of  the  strain 

is  the  energy  used  for  crack  extension  £the  work  of  fracture"] 

W  is  the  energy  released  from  the  system  by  crack  extension 
s 

c  is  a  fracture  length  parameter. 

The  instability  criterion  (Equation  2)  then  becomes 

3(AP)  =  3  [W,  -  W  ]  <  0  (6) 

3c  3c  a  s 

Wd>  the  work  of  fracture,  may  now  include  effects  such  as.  (test) 
strain  rate, tenperature ,  rock  inhomogeneity ,  etc. 

W  may  be  affected  by  variation  of  elastic  modulus  with  stress 
level,  partial  inelasticity,  etc. 

Understanding  of  the  mechanics  of  rock  fracture  and  rational  inter¬ 
pretation  of  rock  strength  tests,  depend  on  the  elucidation  of  the  terins 
'W^  and  W„  for  each  test  situation.  [So  far,  only  the  tensile  fracture 
situation  has  been  considered.  It  will  be  shown  later  that  the  same  remarks, 
apply  for  compressive  strength  tests.] 

It  is  from  a  thorough  study  of  these  terms  that  the  intrinsic  material 


properties  will  be  deduced 
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Stress-Straln  Representation  of  Griffith  Theory  of  'tensile  Rupture 

tee  results  of  laboratory  tests  are  frequently  depicted  as  stress-strain 
curves.  It  is  convenient,  therefore,  to  represent  the  Griffith  theory  in 
tenns  of  the  equivalent  stress-strain  behaviour.  This  has  been  done  for 
the  direct  tension  situation  of  Figure  1,  by  Berry,  [see  publication  (3)  I 
The  result  is  shown  In  Figure  3. 


Figure  Jtress.-s.train  representation  of  Griffith  criterion  of 


ten-. lie  fracture. 


Notes  I'ncracked  plate  1  ui: *  strength  equal  to  theoivtle.a!  col  tec.  1 v«>  trenglh. 
Stress-strain  slope  varies  frcnn(i)  for  uncracked  plate  through  (11),  (111)... 
etc.  for  Increasing  crack  length,  -o-o  indicates  complete  stress-sera  In 
curve  for  slow  extension  of  crack,  to  complete  collapse.  Intial  length  of 
crack  represented  by  slope  (iv). 

As  mentioned  previously,  the  Griffith  strength  depends,  on  the  crack  length 
c,  and  is  not  influenced  by  the  stiffness  of  the  applied  load,  'lliis  mean;:  tint 
the  fracture  initiation,  represented  (for1  cracks  of  various  lengths)  by  the 
locus  ABCD  is  not  affected  by  tlie  load  stlflYiess.  Propagation  of  fracture  Is 
dependent  on  load  stiffness,  'ttius  a(n  unloading)  stl  l’fness  represented  l\v 
the  slope  CQ  will  allow  slow  contrelled  crack  extension  to  specimen  dis¬ 
integration.  A  stiffYiess  as  represented  by  the  slope  OH  will  result  In 
rapid  specimen  disintegration.  Hie  locus  CD  will  not  be  revealed.  Investi¬ 
gation  of  the  physical  processes  of  rock  disintegration  (for  comparison  with 
the  Griffith  model)  requires  that  the  fracture  locus,  be  followed  in  a  control¬ 
led  fashion.  Considerable  attention  has  been  paid  to  this  problem  as  part 
of  Contract  H0101610,  and  a  variety  of  novel  techniques  developed.  Successful 
control  of  the  development  of  fracture  by  the  use  of  soivo-control  testing 
equipment  U:  one  of  the  major  achievements  of  research  under  tills,  contract 
[see  publications  (b),  (8),  ('))]. 

Griffith  Tlieory  of  Rock  Pal  lure  in  Conpressior 

Griffith  attempted  to  extend  the  theory  for  tensile  iupture  to  collapse 
of  specimens  in  compression.  He  assumed  that  fracture  would  occur  when  the 
same  tensile  stress  concentration  (o^)  as  occurred  in  tensile  fracture  was 
reached  at  the  tips  of  cracks  in  specimens  subjected  to  non-liydrostatic  com- 
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press  ive  loading.  The  predicted  corrpressive  strength  was  a  function  of  the 
applied  principal  cotrpressions ,  with  a  predicted  uniaxial  compressive 
strength  8  times  as  great  as  the  tensile  strength.  In  this  extension  of 
the  theory  Griffith  considered  only  the  (necessary)  condition  of  high  tensile 
stress,  ind  (neglected)  the  associated  requirement  of  unstable  energy 
change  with  crack  extension.  Subsequent  research  has  shown  that  tensile 
cracks  do  develop  during  compression  but  that  growth  of  these  cracks  is  a 
stable  process ,  and  the  specimen  is  able  to  support  increased  applied 
loads.  Collapse  of  the  specimen  does  not  occur  until  loads  of  the  order 
of  twice  those  at  which  cracking  begins  are  applied.  The  mechanism  of 
collapse  in  compression  appears  to  be  largely  the  result  of  shear  sub¬ 
sequent  to  the  stable  cracking  of  the  specimen  [see  publication  (1)]. 

Griffith's  neglect  of  the  energy  instability  condition  for  compressive 
loading  thus  resulted  in  an  erroneous  criterion  of  'collapse' . 

Cook  has  developed  a  criterion  of  conpressive  collapse  based  on  an 
■nor,-;/  instability  analysis,  assuming  collapse  in  shear.  This  appears- 
to  b-  an  appropriate  counterpart  to  Griffith's  tensile  fracture  criterion. 

ihe  compressive  strength  (in  plane  stress)  Is  derived  from  the  criterion 


wrier-  r  l.-.  the  .‘.hear  stress  along  the  plane  of  shear’  displacement 
o  is  the  normal  stress  along  the  plane  of  shear  displacement 
u  is  the  coefficient  of  sliding  friction  along  the  shear  plane 
Ys  is  the  "specific  work  of  shear  fracture"  (likely  to  be  affected 
by  axial  splitting) 

c  is  the  original  length  of  the  slrear  crack. 
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Several  points  should  be  noted 

p 

i)  the  (stress)  .(crack  length)  ■  const,  relationship, 

(corresponding  to  Equation  *0  is  again  obtained,  the 

result  of  using  an  elastic  instability  analysis 

ii)  two  quantities  of  energy  supply  and  demand,  similar  to 

W  and  W  , ,  are  involved 
s  d* 

ill)  the  criterion  may  be  presented  as  a  stress-strain  locus 
similar  to  Figure  3. 

'Hius,  man/  of  the  concepts  used  in  tensile  fracture  may  also  be  applicable 
to  compressive  collapse. 

Experimental  results  reveal  that  controlled  collapse  of  rock  specimens 
in  compression  does  result  in  a  complete  stress-strain  curve  of  the  general 
form  predicted  in  Figure  3,  although  the  development  of  macroscopic  shear 
■ through  interaction  with  cracking)  is  more  complicated  than  assumed  in 
Cook's  model,  [see  publication  (10)3 


Size-Strength  Relationship  -  (predictions  of  Griffith  theory) 

Hie  form  of  tire  energy  change  equation.  Equation  (1)  on  which  Griffith's 
tiieory  is  basod,  has  an  interesting  implication  concerning  the  effect  of 
change  in  specimen  size  on  the  specimen  strength.  Equation  (1)  is  expressed 
3*3 


or 


AP  *  ^cy  - 
AP  *  Wd  - 


2  2 
TtC  O 


(1) 


Note  that  the  first  tem,  the  ' energy  demand*  tern  (W^)  has  the  dimensions 
of  area  (since  the  plate  thickness  is  a  unit  length)  whereas  the  second. 
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'energy  supply'  tenn  (W  )  has  the  dimensions  of  volume.  If  the  crack 

s 

dimensions  increase  then  the  strength  will  decrease,  simply  because  the 

rate  of  increase  of  W  with  size  exceeds  that  of  W ..  This  is  in  fact 

s  d 

wtty ,  from  the  energy  stability  analysis,  a  dependence  of  strength  on  size 

of  crack  such  as  indicated  by  Equation  4,  is  obtained. 

Although  Equation  (1)  strictly  applies  to  an  Infinitely  large  plate, 

the  d  imens  ional  difference  between  W.  and  W  will  reran  In  even  for  finite 

Cl  s 

test  specimens,  and  for  Interacting  cracks. 

Tliis  has  several  interesting  implications  concerning  the  stability 
of  large  structures  as  well  as  small  structure,  and  lias  been  discussed 
in  publication  (2).  Some  inferences  into  the  possible  existence  of  a 
size-strength  dependency  can  be  drawn  from  examination  of  the  dependence 
of  ’c'  in  Equation  (1)  or  Equation  CO  on  specimen  size. 

Hie  length  'c'  characterizes  tlie  most  severe  flaw  (considering 
level  of  stress  intensity,  length,  and  orientation)  in  the  specimen. 

The  remaining  flaws  combine  with  the  'intact'  portions  of  the  specimen 

to  give  the  overall  modulus  E  of  the  'plate'  in  Figure  1.  If  a  larger 

specimen  (plate)  is  chosen,  then  there  is  some  probability  that  it  will 
contain  a  larger  more  severe  'crack'.  In  the  extreme  case  of  an  extensively 

cracked  material,  the  maximum  crack  size  in  a  specimen  may  be  almost 

directly  pronortional  to  the  size  of  the  (geometrically  similar)  specimen. 

If  we  characterize  the  specimen  by  a  length  L,  then  we  have 

c  a  E  (8) 

Equation  (H)  may  then  be  rewritten  as 

4-  r  -0.5 
o  =  const.  L 
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A  relationship  of  this  form  has  been  found  to  hold  approximately  for 

the  strength  of  geometrically  similar  coal  pillars.  Coal  probably 

comes  closer  than  any  other  'rock'  to  satisfying  Equation  (8),  being 

characterized  by  extensive  internal  cracking. 

If,  on  the  other  hand,  we  assume  that  a  specimen  is  composed  of 

a  large  number  of  identical  flaws  (requiring  the  same  Wd)  equidistant 

from  each  other,  with  an  elastic  region  around  each  flaw  supplying 

the  energy  (W  ) ,  then  the  strength  of  each  flaw  will  be  the  same . 
s 

Increase  in  specimen  volume  will  not  affect  the  strength,  i.e.  Equation 
(4)  becomes 

o  *  const.  L°  (9) 

Ihe  'flaws'  in  this  case  may  be  related  to  the  grain  siz  in  a 
homogeneous  material  [such  as  the  Georgia  Cherokee  marble  used  in  size 
and  shape  effect  studies.  See  publication  (7)j 

Thus,  one  may  make  an  estimate  of  the  size-strength  dependency  from 
the  degree  to  which  large  'flaws'  are  introduced  in  larger  specimens. 

^he  exponent  'o'  in  the  strength-size  (constant  shape)  relationship 

o  *  const.  La  (10) 

may  vary  over  the  range  0  <  a  <  -0.5  with  the  lower  extreme  (a  =  -0.5) 
being  approached  for  extensively  cracked  material. 

Strength  of  Jointed  Materials 

The  above  discussion  suggests  that,  if  a  specimen  is  composed  of  joints  of 
essentially  Identical  properties  (Wd  =  const.),  then  the  strength  will 
depend  on  the  intact  volume  associated  with  supply  of  energy  (Wo)  to  each  joint. 
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Ln  other  wonis  the  strength  will  depend,  within  limits,  on  tiie  spacing:  between 
Joint:;.  At  the  lower*  limit,  of  very  closely  spaced  Joint.:;,  the  energy 
supply  will  decline  proportionately  due  to  interaction  between  adjacent 
joints.  At  the  upper  limit,  widely  spaced  joints,  the  energy  supply  (W  ) 

(i.e.  as  released  with  increment  joint  shear)  will  tend  towards  a  finite 
value.  The  strength  of  a  jointed  specimen  should  thus  exhibit  a  dependence  on 
joint  spacing  over  a  range  of  spacing.  Experiments  to  test  this  belief  were 
planned  for  the  early  months  of  the  second  year  of  Contract  H0101610. 

Size  Effect  in  Ring  Tests 

The  above  discussion  has  implicity  assumed  that  the  smallest  specimen 
is  large  relative  to  the  grain  size,  and  that  the  specimen  is  homogeneously 
stressed.  In  some  instances  (e.g.  in  indirect  tension  tests)  the  critically 
stressed  region  is  a  small  part  of  the  specimen  and  the  applied  stresses 
may  vary  considerably  over  the  length  of  a  grain.  In  such  cases,  although 
the  energy  instability  mechanism  is  probably  valid,  the  boundary  conditions 
of  the  Griffith  plate  model  are  far  from  satisfied  and  deviations  in  strength 
may  be  expected.  The  problems  of  tensile  strength  determination  using  indirect 
methods  have  been  extensively  studied  in  Publication  (*J).  In  sunmary  it 
may  be  said  that  most  such  tensile  tests  are  conducted  with  test  specimens 
that  are  much  too  small.  Results  are  currently  being  re -examined 
from  the  energy-instability  viewpoint. 

Energy  Instability  Analysis  of  Indirect  Tests 

Analytical  computation  of  the  energy  changes  associated  with  crack 
extension  in  an  applied  stress  field  is  extremely  difficult  except  for 
trie  very  simple  conditions,  assumea  by  Griffith.  The  introduction  of  the 


Finite  Element  Modal  of  Olotnotrol  Comprmion 
Ring  Toal,  With  Vortical  Crock 


FINITE  ELEMENT  REPRESENTATION  OF  A  DIAMETRALLY 
LOADED  DISC  WITH  A  SMALL  CENTRAL  HOLE. 


VCOTCM.  ONPUKtMtNT  I, 

Homwi^M.  ooruuncNT 


"GRIFFITH  LOCUS"  FOR  A  DIAMETRALLY  LOADED  DISC 
CONTAINING  A  SMALL  HOLE. 


Figure  4 .  Griffith  locus  for  a  dlametrally  loaded  disc  containing 


a  small  hole. 
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finlte  element  method  has  facilitated  the  study  of  energy  change  - 
crack  extension  relationships  for  finite  size  of  specimens  in  complex 
stress  fields. 

A  notable  accomplishment  under  Contract  H0101610  has  been  the 
successful  application  of  the  finite  element  method  to  determine 
the  energy  changes  with  crack  extension  for  indirect  tensile  tests. 

The  method  is  outlined  in  Publication  (3).  The  complete  force-displacement 
(vertical  and  horizontal)  curve  for  the  ring  test  has  been  derived,  and 
is  shewn  in  Figure  *1  compares  well  with  experimentally  observed  curves. 

The  method  is  applicable  to  other  tests  and  is  currently  being  used  to 
analyze  beam  test  data.  The  results  will  be  reported  in  Publication  (12) 
now  in  preparation. 

It  is  confidently  expected  that  the  study  of  -indirect  tensile  tests 
in  this  way  will  clear  up  much  of  the  confusion  that  now  exists  concerning 
the  tensile  strength  of  rock  specimens.  It  should  also  establish  the 
nature  of  the  intrinsic  nulerial  properties  which  determine  the  tensile 
strength  of  rocks. 

Influence  of  End  Restraint  In  Compression  Testing 

Rock  compression  tests  are  usually  interpreted  on  the  assumption  that 
the  load  is  uniformly  distributed  across  the  specimen  at  all  stages  of  the 
test.  It  Is  v/ell  recognized  that  this  can  not  be  true  since  the  properties 
of  the  platen::  remain  elastic  whereas  those  of  the  specimen  change  throughout 
the  loading  cycle.  Hie  extent  to  which  end  restraint  influences  compressive 
failure  has  been  studied  in  four  ways  under  Contract  H0101610  as  follows 


1.  A  special  sample  configuration  (viz.  hollow  cylinders  with 
geometrically  matched  platens)  which  substantially  reduces  end  con¬ 
straint  (demonstrated  from  a  finite  element  study)  has  been  used  to 
test  rock  cylinders.  Hie  specimens  were  sectioned  and  the  develop¬ 
ment  of  fracture  examined.  It  was  concluded  that  end  restraint  does 
not  significantly  change  the  indicated  compressive  strength  or  the 
compressive  modulus.  Hie  findings  of  this  study  sire  reviewed  in  detail 
in  Publication  (1). 

2.  A  theoretical  and  experimental  study  of  the  influence  of  lateral 
stiffness  of  the  machine  on  specimen  collapse  in  compression  reveals 
that  this  (previously  unconsidered)  stiffness  may  exert  significant 
restraints  on  macroscopic  shear  of  specimens,  and  could  thereby  affect 
compressive  strength.  A  laterally  stiff  system  also  tends  to  Indicate 

a  pseudo  size  strength  effect  where  none  actually  exists.  This  effect 
will  be  fully  discussed  in  the  forthcoming  publication  (13).  Methods 
of  testing  to  determine  the  influence  of  lateral  stiffness  will  be 
described.  (See  Appendix  2  for  abstract) 

3.  Hollow  cylinders  (5  inch.  O.D.  1  inch  I.D. )  of  rock  were  subjected 
to  increasing  hydrostatic  external  pressure  until  compressive  collapse 
occurred  on  the  inner  diameter.  Final  collapse  usually  occurred  in  macro¬ 
scopic  shear.  Procedures  and  results  are  reported  in  E’ublication  (2). 
Indicated  compressive  strengths  for  this  mode  of  loading  were  200%  -  300% 
greater  than  the  values  obtained  by  the  standard  axial  compression  of 
cylinders.  This  indicates  that  the  compression  'strength',  like  tensile 
'strength'  is  not  a  true  material  constant.  It  seems  very  likely  that 
the  variation  is  due  to  the  difference  in  stress  gradients  (i.o.  energy 
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density)  developed  in  the  two  tests.  Energy  instability  analysis 
should  help  clarify  the  cause  of  the  discrepancy. 

Further  tests  with  a  larger  O.D.  pressure  cylinder  and  specimens 
are  recorrmended  to  better  avoid  grain  size  inhomogeneity  effects.  The 
pressure  vessel  should  be  designed  to  withstand  50,000  -  60,000  p.s.i. 
in  order  to  cause  compressive  collapse  of  the  specimens. 

An  instrumented  loading  head  has  been  designed ,  constructed  ;md 
instrumented  to  enable  the  normal  load  distribution  at  the  rock-platen 
interface  to  be  monitored  throughout  an  axial  compression  test.  A 
finite  element  program  has  been  written  to  enable  the  interface  stresses 
to  be  deduced  from  the  stresses  measured  by  the  loading  head.  This 
instrument  will,  for  the  first  time,  enable  the  applied  load  distribution 
to  be  examined  in  detail  and  correlated  with  fracture  development. 

Questions  such  as  those  raised  in  Publication  (7)  concerning  the  actual 
load  bearing  area  of  a  conpressive  specimen  under  load,  and  correction 
to  the  'true  stress-strain'  curve  should  be  answered  by  the  information 
provided  by  this  loading  head. 

Experimental  Determination  of  Complete  Stress-Strain  Curves  in  Compression 
and  Indirect  Tens  ion 

The  successful  development  of  experimental  methods,  for  obtaining  complete 
stress-strain  curves  in  compression  and  indirect  tension  tests  even  for  tin: 
case  of  Class  II  failure  [see  lublication  ( 1 1 )  J  is  a  very  significant  step 
towards  the  goal  of  determining  the  intrinsic  material  properties.  Further 
attention  must  be  given  to  evaluating  the  influence  of  specimen  test 
conditions  on  the  observed  complete  'stress-strain'  curve  before  the  latter 


STRESS  p •  • 
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can  be  used  in  the  derivation  of  real  material  properties.  Complete  details 
of  progress  to  date  are  given  in  Publication  (9).  The  observation  of  the 
effect  of  strain  rate  on  the  complete  stress -strain  curve  could  iiave 
important  Implications  concerning  the  relatively  unique  nature  of  th 
post-peak  load  curve.  Figure  5  discussed  in  Publication  (5)  reveals 


Figure  5 .  Time  deformation  at  constant  displacement. 

that,  even  though  deformation  was  periodically  halted  for  one  hour  (at  each 
"Hold")  the  load-deformation  path  followed  on  resumption  of  deformation 
[strain  rate  1.5  x  10”^/second]  was  essentially  the  same  as  that  followed 
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during  continuous  defomatlon.  Other  related  studies  also  suggest 
that  the  post-peak  curve  may  have  some  fundamental  significance  in 
creep  testing.  It  may,  for  example,  mark  the  strain  in  a  creep  test 
(at  a  certain  stress  level)  at  which  tertian  creep  ensues.  Such  a 
finding  would  liave  considerable  implications  in  tiie  design  of  under¬ 
ground  mine  pillars. 


Lefomablllty  of  Rock  Specimens 
Modulus  of  Deformation 

The  defomability  of  (Isotropic)  rock  specimens  is  usually  expressed 
in  terns  of  a  modulus  of  deformation  (elasticity)  and  Poisson's  ratio. 
Experimentally,  the  modulus  of  deformation  is  measured  as  the  ratio  of  the 
average  load  per  unit  cross-sectional  area  of  the  specimen  to  the  overall 
contraction  (between  platens)  per  unit  length  of  the  specimen.  [Note: 
strain  gauges  attached  to  detect  axial  strains  over  the  central  part  of  a 
specimen  give  an  erroneous  measure  (of  around  7%)  of  the  modulus  if  related 
to  the  mean  axial  compression.]  The  effect  of  platen  end  restraint  increases 
as.  tiie  height  to  width  ratio  of  the  specimen  is  decreased.  A  simple  cal¬ 
culation  of  two  extreme  cases  where  (a)  no  end  restraint  exists  and  (b) 
where  eni  restraint  is  such  as  to  completely  inhibit  axial  expansion  [l.e. 
e  =  r  =  0  (if  2  i:;  taken  as  the  axial  direction,  nils  condition  is 

v  •  t 


approached  for  very  squat  specimens.]  indicates  that  the  ratio  of  the 

apparent  modulus  of  elasticity (E,  )  to  the  tree  modulus  (E)  varies  iVoni 

app 

E  =  1.0E  for  case  (a),  to 

app 

E  -  (l~v)  E  for  case  (b) 

aPP  (l+v)Tl-2v7 
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Thus  if  v  is  taken  as  0.25  for  most  rocks  over  the  elastic  range  of  deformation 

E  =  1.P0E  for  case  (b) 
app 

Thus,  a  maximum  error  of  +20%  may  occur  in  the  estimation  of  K  from  the 
ratio  of  mean  stress  to  axial  deformation.  In  most  real  situation:-  end 
constraint  is  less  severe  and  the  error  will  be  reduced.  A  finite  element 
analysis  of  the  effect  of  end  restraint  on  the  error  in  determination  of 
the  modulus  of  elasticity  is  currently  being  made.  It  seems  very  probable 
that  results  will  confirm  the  experimental  finding  [see  Publication  (7)] 
that  the  observed  modulus  of  elasticity  is  virtually  uninfluenced  by  end 
restraint. 

This  pen, Its  the  important  conclusion  to  be  drawn  that  the  modulus 
of  elasticity  measured  in  conpression  can  be  considered  to  be  a  true 
material  property.  It  should  be  noted k.  however,  that  the  modulus  may 
change  with  level  of  loading  and  also  will  probably  differ  for  loading 
and  unloading  conditions.  A  similar  conclusion  may  be  drawn  concerning 
the  modulus  in  tension.  However,  it  should  be  noted  that  the  modulus  in 
tension  will  generally  not  be  the  same  as  the  modulus  in  compression  and  should 
be  determined  directly. 

b)  Poisson's  ratio 

No  work  has  been  conducted  as  yet  on  the  value  of  To  is, son's,  ratio.  The 
uiyeement  between  the  Poisson's  ratio  of  the  material  and  that  comput'  d  iVom 
inoasurement  of  the  lateral  ;ind  axial  displacement  will  depend  on  the  degive  of 
internal  fracturing  within  the  specimen.  Pre-existing  open  axial  fractures 
will  result  in  under-estimation  of  Poisson's  ratio.  Axial  cracks  developed 
during  loading  will  result  in  over-estimation  particularly  at  loads  approach¬ 
ing  the  maximum  and  in  the  post-peak  loading  range,  where  the  (elastic)  notion 
of  Poisson's  ratio  is  entirely  invalid. 
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DEF0RMAB1UTY  AND  STRfcW&TH  OF  ROCK  MASSES 

The  ultimate  aim  of  the  research  started  under  Contract  H0101610 
is  to  derive  rational  methods  for  estimating  the  deformability  and  strength 
of  rock  masses  subjected  to  load  as  a  result  of  engineering  activities. 

Emphasis  in  the  first  year  has  been  placed  on  a  study  of  the 
deformability  of  rock  masses  rather  than  on  strength.  The  main  reason  for 
this  is  the  one  already  rioted  (i.e.  the  basic  concept  that  the  strength 
defines  the  onset  of  unstable  energy  changes  should  be  as  equally  valid  for 
reek  masses  as  for  rock  specimens),  it  is  preferable  to  delay  Intensive 
stud;,  of  the  applicability  of  these  ideas  to  rock  mas, res  until  they  have 
been  fully  explored  and  tested  on  laboratory  specimens.. 

DEFORMABIL  LTY  OF  ROCK  MASSES 

'Hie  rational  understanding  of  the  deformabilitv  of  a  rock  mass  is 
Important  in  many  engineering  problems.  These  Include,  for  example, 
the  problem  of  dam  design,  where  the  deformability  of  the  lock  foundation 
is,  important,  water1  pres, sure  tunnels  where  the  deformability  of  the 
tunnel  wall  rock  influences  lining  design,  and  surface  construction  where 
the  deformability  of  the  foundation  may  affect  the  stresses  developed  In 
tt.  >  surface  structure. 

Studies  of  the  deformability  of  rock  masses  under  Contract  HOlOltdO 
luncheon  restricted  entirely  to  theoretical  analyses  of  the  general  problem. 

Early  in  the  stud;/  visits  were  made  to  industrial  research  laboratories 
in  an  attempt  to  obtain  atcual  field  data.  Although  some  valuable  insight  wa 


obtained  into  current  design  of  engineering  structures  of  which  mass  deform- 
ability  is  important ,  it  was  quickly  realized  that  much  of  the  needed  data 
was  simply  not  available.  'ihis  is  due  in  part  to  the  absence  of  a  theoretical 
framework  from  which  to  plan  side  investigations.  It  was  therefore  decided  to 
develop  a  finite  element  model  to  simulate  a  jointed  rock  mass  on  which 
simulated  tests  could  then  be  performed.  The  model  w;is  mite  general  and 
allowed  a  wide  variety  of  practical  rock  mass  situations  to  be  simulated. 

In  this  way  field  data  could  synthetically  be  generated  for  use  in  cheeking, 
proposed  methods  of  analyzing  (real)  field  results,  'ihe  method  has  the 
advantage  that  the  full  scale  properties  of  the  mass  were  known,  so  that  the 
accuracy  of  any  measurement  technique  could  be  readily  be  checked. 

’Ihe  analytical  work  to  date  has  concentrated  on  the  development  of 
rational  programs  for  determining  in-situ  rock  mass  defonnabillty .  'ihe 
main  experimental  methods  considered  were : 

1.  Radial  jacking  tests 

2.  Borehole  dilatometer  tests 

3.  Plate  loading  tests. 

A  variety  of  serious  theoretical  obstacles  lias  to  be  overcane  in 
progressing  towards  this  practical  goal.  Considerable  effort  wa-  extended, 
for  example,  in  obtaining  an  analytical  solution  for  the  mean  displacement 
of  a  pressurized  dilatometer  when  situated  in  a  cylindrical  borehole 
intersecting  a  planar  jointed  rock  mass.  This  is  a  three  dimensional  situation 
which  is  completely  impractical  to  solve  by  finite  element  analysis  using 
currently  available  cauputers.  'Ihe  problem  was  solved  and  narks  a  very 
significant  achievement  of  the  research.  This  and  other  theoretical  results 
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wiil  tie  discussed  in  detail  a:?  pai’t  of  future  publications,  some  of  which 
are  outlined  in  the  abstracts  appended  to  this  report.  Some  details  have 
already  been  given  in  the  quarterly  progress  reports  submitted  as  part  of 
Contract  H0101610. 

For  the  purposes  of  this  report  the  year's  progress  towards  the 
practical  goal  of  providing  guidance  in  the  rational  planning  and  inter¬ 
pretation  of  deformability  measurements  will  be  chronologically  outlined. 
Following  niis  however  the  reader  is  reminded  that  many  of  the  analytical 
achievements  of  this  work  are  glossed  over. 

Significance  of  Jointed  Character  of  Rock  Masses 

'Hie  first  question  to  be  resolved  in  order  to  characterize  a  jointed 
rock  mass  is  the  following: 

How  important  is  it  to  accurately  represent  the  detailed  jointed 
structure  of  a  rock  mass  in  order  to  estimate  the  mass  properties? 

Is.  it  sufficient  for  example  to  make  the  assumption  [as  war,  done 
by  Cu  et  al.  (Eleventh  Symposium  on  Rock  Mechanics,  253-266)]  that  a  rock 
mass  could  be  represented  as  a  finite  element  network  in  which  the  isotropic 
elastic  properties  of  each  element  vary  randomly?  Finite  element  analysis 
curled  out  ho  examine  tills  question.  'Use  rock  muss  was,  divided  Into 
unit  (elemental)  masses  and  the  isotropic  elastic  properties  of  each  unit, 
allowed  to  vary  randomly,  'die  finite  element  grid  was  then  super  imposed 
on  the  truss.  The  properties,  of  eacii  (variable  size)  element  of  ttie  Unite 
element  >7dd  were  taken  as  trie  mean  value  of  tiie  properties  of  the  unit 
contained  in  the  element.  Hie  influence  of  the  shape  of  the  unit  on  the 
t  repo  sties  of  the  finite  element  was  also  Included.  Thus  the  mass  was 
con:'  1  Jered  to  be  represented  in  each  of  three  ways : 
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a)  As  a  series  of  semi-infinite  layer  units  (on  a  rigid  base  at  depth) 

b)  As  a  series  of  cylindrical  units 

c)  As  a  series  of  spherical  units 

Three  independent  finite  element  programs  were  developed  to  determine  the 
deflection  of  a  point  load  applied  normally  to  the  surface  of  an  elastic 
mass  (i.e.  the  Boussinesq  problem)  composed  of  units  as  described  for  the 
tliree  cases  above.  Each  situation  was  analyzed  as  a  randan  heterogeneous 
inclusion  problem.  The  properties  of  each  of  the  units  were  randomly 
selected.  The  tests  were  then  repeated  for  a  different  set  of  random  pro¬ 
perties  (i.e.  the  Monte  Carlo  simulation  technique  was  used)  and  the 
point  displacements  again  determined.  The  coefficient  of  variation  (i.e. 
the  ratio  of  standard  deviation  to  the  mean  displacement)  was  then  determined 
for  each  of  the  rock  mass  models  a,  b  and  c  above.  The  following  results 
were  obtained. 

a)  semi-infinite  layer  model,  coefficient  of  variation  1  to  0.1, 

b)  cylindrical  rod  model,  coefficient  of  variation  0.1  to  0.01, 

c)  spherical  model,  coefficient  of  variation  0.01  to  0.001. 

'Ihus,  the  coefficient  of  variation  ranged  over  three  orders  of  magnitude 
depending  on  the  model  assumed  for  the  elements  of  the  mass. 

Tills  study  demonstrated  conclusively  that  it  is  important  to  correctly 
model  the  geometry  in  finite  element  analyses  of  rock  masses  if  realistic 
results  are  to  be  obtained,  [see  2nd  Quarterly  He  port  | 
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The  jointed  rock  mass  was  also  represented  as  a  semi -infinite  rock  mass 
containing  randoinly  oriented  joints  of  random  normal  and  shear  stiffYiesses. 

Only  two  dimensional  joint  sets  were  considered.  The  joint  planes  were 
also  normas  to  the  section  through  the  joint,  and  the  coefficient  of  variation 
in  spacing  orientation  and  stiffness  did  not  exceed  the  wide  limit  of  50% . 

Test  data  were  then  obtained  for  plate  bearing,  dilatometer  and  radial 
jacking  tests.  By  computing  the  mean  displacements  normal  to  the  applied 
test  load  the  instrument  displacement  was  computed  by  determining  the  increase 
in  strain  energy  of  the  mass  due  to  the  introduction  of  the  joint. 

Sets  of  displacement  values  were  obtained  by 

a)  varying  the  assumed  joint  system, 

b)  varying  the  position  of  point  of  load  application  with  respect 
to  the  joint, 

c)  varying  the  size  of  the  loaded  area. 

'Ll i is  enabled  tire  distribution  of  measured  displacements  to  be  determined 
for  each  size  of  instrument.  From  this  the  mean  standard  deviation  and 
confident  intervals  were  computed  for  a  fixed  size  of  instrument  and  a  given 
rock  trass.  The  number  of  tests  required  to  obtain  a  mean  displacement  or 
equivalently  modulus  of  deformation  within  plus  or  minus  30 %  of  the  true  mean  mass 
displacement  at  <)0%  confidence  level  was  then  determined.  Tire  results  of  all 
the  rock  analyses  can  be  combined  into  a  single  curve.  (Bee  Figure  6). 

Figure  6  implies  that  the  difference  between  the  modulus  ;md  deformation 
measured  by  radial  jack  tests  and  dilatation  tests,  will  not  differ  by  more 
than  30?,  provided  that  the  orientation  of  both  Instruments,  is  the  same,  and 
the  number  of  tests  conducted  is  as  demanded  by  Figure  6  shown  on  the  next  page. 
Obviously,  since  the  pressurized  region  in  a  dilatometer  is  usually  much 
shorter  than  in  a  plate  jacking  test,  the  required  number  of  dilatometer  tests 


NUMBER  OF  TESTS 
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AVERAGE  NUMBER  OF  JOINTS  INTERSECTING  THE  DILATOMETER 
[OR  RADIALLY  (JACK)  PRESSURIZED  REGION]  PER  TEST 


Figure  6.  Number  of  dllatometer  tests  required  to  estimate 
a  modulus  of  reduction  factor  in  a  jointed  rock 


mass. 


will  be  considerably  higher. 

The  development  of  the  curves  in  Figure  6  is  considered  to  represent 
a  very  significant  step  towards  the  stated  goal  of  establishing  rational 
procedures  for  the  determination  of  tire  deformation  modulus  of  a  jointed 
rock  mass.  It  is  likely  that  the  model  will  need  to  include  practical 
detail::,  not  yet  considered  but,  for  the  first  time,  a  method  is  presented 
which  will  suggest  a  number  and  orientation  of  tests  necessary  to  determine 
the  mass  modulus  to  a  specified  degree  of  accuracy. 
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It  is  particularly  interesting  that  the  results  appear  to  suggest  that 
the  number  of  dilatometer  tests  required  to  determine  the  modulus  with 
an  accuracy  comparable  to  that  of  a  single  plate  jack  test  is  usually 
such  that  the  dilatometer  method  would  be  economical. 

'Hie  abstract  of  the  three  proposed  papers  by  Dr.  Singh  includes 
several  specific  findings  related  to  the  above  mentioned  general  results. 

In  summary,  substantial  progress  has  been  made  toward  the  objectives 
of  defining  necessary  test  programs  to  ensure  accurate  measurements  of 
the  defonnability  of  rock  masses.  Versatile  analytical  procedures  are 
fully  operational.  They  can  readily  be  used  to  study  laboratory  specimens 
of  rock.  Close  attention  can  be  given  to  logical  situations  which  are 
likely  to  arise  and  can  be  adequately  represented  by  this  theoretical 


model . 
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LliVl’  OF  MAJOR  ACHIEVHKMTS  IN  ARPA  CONTRACT  110101610 

1.  The  basic  concept  of  the  Griffith  theory  (i.e.  for  failure  to 
occur  there  must  be  a  balance  between  the  rate  of  energy  required 
and  the  rate  of  energy  available)  has  been  extended  to  a  variety 
of  situations  including  laboratory  tests,  tunnel  collapse  and 
blasting.  The  concept  is  also  equally  applicable  to  the  analysis 
of  rock  mass  stability  and  strength. 

2.  Theoretical  models  (incorporating  finite  element  techniques)  have 
been  developed  for  constructing  the  complete  Griffith  locus  and 
complete  force-displacement  curve  for  any  type  of  test  on  rock  in 
compression  or  tension. 

3.  Completely  new  servo-control  techniques  have  been  developed  for  experi¬ 
mentally  obtaining  the  complete  force-displacement  curve  in  any  type  of 
rock  mechanics  test.  These  have  enabled  detailed  studies  of  rock  collapse 
in  many  laboratory  tests. 

4.  From  both  theoretical  and  experimental  results  it  can  be  definitely 
stated  that  the  strengths  of  rock  (tensile  and  compressive)  are  not 
material  properties  and  that  the  strength  and  failure  behavior  of 

a  rock  specimen  or  rock  mass  can  only  be  derived  using  the  minimum 
potential,  energy  type  of  approach. 

'j.  'Die  influence  of  end  restraint  during  compressive  testing  has  been 

studied  in  four  ways  and  the  design,  construction  and  instrumentation 
of  a  loading  head  that  enables  the  actual  stress  distribution  applied 
to  the  specimen  ends  to  be  determined  at  all  stages  of  the  complete 
stress-strain  curve  represents  a  major  advance  towards  complete 
landerstanding  of  end  restraint. 
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otudies  or  the  time  dependent  behavior  of  failed  rock  indicate  tliat 
the  complete  stress-strain  curve  is  strongly  time  dependent  (more 
so  than  intact  rock  because  the  specimen  is  in  a  continual  state  of 
progressive  collapse)  and  that  the  time  factor  must  be  included  in 
any  analysis,  of  the  strength  and  failure  behavior  of  rock  masses. 

It  lias  been  conclusively  established  that  the  measured  initial  elastic 
modulus  is  a  material  property.  Sigiif  leant  variations  could  occur 
in  unusual  geometrical  testing  situations  (very  low  specimen  height: 
width  ratios)  but  for  normal  laboratory  testing  the  measured  modulus 
should  only  vary  from  the  true  modulus  by  a  maximum  of  a  few  per 
cent . 

A  finite  element  model  has  been  developed  in  quite  general  form 
that  allows  a  wide  variety  of  practical  rock  mass  situations  to 
be  simulated. 

An  analytical  solution  has  been  obtained  for  the  mean  displacement 
of  a  pressurized  dilatcmeter  when  situated  in  a  cylindrical  borehole 
intersecting  a  planar  Jointed  rock  mass  in  theee  dimension.  This 
marks  a  very  significant  achievement  of  the  research. 

Finite  element  analyses  of  rock  masses  were  made  in  which  the  mass 
was  divided  into  different  types  of  elements  and  random  properties 
assigned  to  each  element.  It  was  found  that  the  coefficients  of 
variation  for  the  displacements  in  the  Boussinesq  problem  varied  through 
ttiree  orders  of  magnitude  depending  on  the  assumed  model.  This 
proves-  conclusively  chat  the  geometry  of  rock  masses  must  be 
realistically  simulated  in  finite  element  analyses. 
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11.  'Hio  distribution  ol‘  displacement:-,  was  calculated  for  plate  hearing, 
dllatometer  and  radial  Jacking  tests  for  various  Jointed  system:-, 
and  size  of  the  loaded  area,  liie  number  of  tests  required  to  obtain 
a  modulus  of  deformation  within  certain  confidence  limits  was  then 
determined  and  this  represents  a  high  plateau  of  achievement  towards 
the  stated  goal  cf  establishing  rational  procedures  for  the  determination 
of  the  deformation  modulus  of  a  jointed  rock  mass. 
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Attract  of  Publication  12  -  '"Die  failure  of  rock  beams11  by  J.  A.  HUDSON, 

M.  P.  HARDY  and  C.  FAIRHURST 

It  is  noted  that  there  is  considerable  variation  in  the  tensile  strength 
of  rock  as  determined  by  the  beam  test  (the  modulus  of  rupture).  This 
variation  is  of  three  basic  types:  with  repeated  tests,  with  beam  sise 
and  shape  and  with  other  indirect  tensile  tests. 

Weibull’s  statistical  theory  of  the  strength  of  materials  is  discussed 
:ind  evidence  presented  to  Indicate  Weibull ' s  theory  does  not  adequately 
characterise  the  failure  of  rock.  Complete  force-displacement  curves 
representing  complete  beam  collapse  are  derived  or.  more  mechanistic  grounds 
using  the  basic  Griffith  concept  and  finite  element  analysis.  The  implications 
of  this  analysis  are  discussed  and  compared  with  experimental  data  obtained 
from  tests  in  which  the  complete  beam  collapse  process  was  controlled  and 
photographed  in  a  servo-controlled  testing  machine. 

Ai  struct  of  Publication  13  -  "Influence  of  the  lateral  stiffness  of  the 
testing  machine  -  specimen  system  on  failure  characteristics  of  rocks” 
by  F.  COHHhT  and  C.  R.  NELSON 

A  theoretical  model  Is.  developed  to  describe  the  complete  l’orce- 
displ-i cement  behavior'  of  a  rock  specimen  tested  in  conpresslon  and  failure 
by  Localized  shear.  'Hie  influence  of  the  lateral  stiffness  of  the  testing 
machine  -  specimen  system  on  the  failure  of  specimens  is  analysed. 

For  certain  loading  situations,  the  lateral  stiffness  of  the  machine 
i.  predicted  to  have  an  Influence  on  the  force-displacement  curve  and  the 
tyjx  of  failure  mode:  this  is  in  addition  to  the  specimen  size  and  shape  e'Toct. 


In  particular,  the  analysis  predicts  effects  compatible  with  experimental 
observations  and  the  change  from  asynmetric  shear  failure  to  symnetrical 
shear  failure  for  increasing  lateral  machine  stiffness. 

Properties  determined  in  a  testing  machine  of  zero  lateral  stiffness 
will  be  less  experimentally  sensitive  than  those  obtained  by  current 
methods.  To  overcome  the  difficulties  of  using  such  a  testing  machine,  a 
method  is  proposed  to  determine  the  basic  compressive  behavior  from  the 
properties  of  peak  cohesion,  the  deformation  required  to  destroy  all 
cohesion  and  the  coefficient  of  internal  friction,  Hiese  properties  are 
determined  in  separate  tests  and  can  be  considered  to  be  intrinsic  material 
properties. 

Abstract  of  Put’ll  cation  14  -  "Continuum  characterization  of  a  rock  mass  - 
Part  I:  The  constitutive  equations”  by  B.  SINGH 

Hill  (1963)  has  proved  that  the  average  strain  ener>y  density  in  any 
region  of  an  elastic  and  inhomogeneous  material  can  bo  calculated  from 
the  averages  of  the  stresses  and  strains  within  that  region.  This  concept 
has  been  used  to  derive  the  general  constitutive  equations  of  a  rock  mass 
containing  an  orthogonal  set  of  discontinuous  Joints  intersecting  an 
anisotropic  rock  material.  Hie  constants  required  for  the  continuum 
characterization  of  the  jointed  mass  are  tlie  joint  stress,  concent  rat  Ion 
factors  Bjh  and  (see  Figure  A.l).  'Ihese  are  defined  as  the  ratio  of  sire 

along  the  joint  to  the  overall  stresses  in  the  rock.  In  the  case  of 
continuous  joints,  these  stresses  concentration  factors  are  equal  to  unity. 
Tais  leads  to  simple  anisotropic  stress-strain  relations. 
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FIG  A.l  ROCK  MASS  CONTAINING  STAGGERED  JOINTS 


It  is  also  shown  that  the  anisotropic  stress-strain  relations  for  the 
rock  mass  can  be  computed  using  average  joint  compliances  and  average 
joint  frequency,  when  these  parameters  vary  randomly,  i.e.  it  is  not  necessary 
to  consider  the  entire  range  of  individual  joint  properties.  Further,  exact 
expressions  have  been  obtained  for  the  stress  concentration  factors  in  case 
of  a  rigid  rock  containing  staggered  compliant  joints  and  a  correction  is 
proposed  for  elastic  rocks.  These  results  suggest  that  interlocking  between 
blocks  of  the  rock  may  become  significant  even  for  a  slight  discontinuity 
(i.e.  offset)  along  the  joints  such  as  may  occur  during  shear  deformation 
of  a  rock  mass  with  an  initially  continuous  orthogonal  joint  set.  Stress 
concentration  factors  computed  independently  from  the  results  of  a  finite 
element  program  for  a  jointed  mass  compare  excellently  with  the  above 
mentioned  theoretical  results. 

It  is  further  shown  that  tensile  stresses  are  developed  inside  a  rock 
with  staggered  joints,  and  may  be  as  high  as  twice  the  overall  shear  stresses 
or  the  overall  corrpressive  stresses.  Typical  stress  distributions  within  a 
block  are  shown  in  Figure  A. 2a  and  A. 2b. 

It  is  concluded  that  a  rock  mass  is  rendered  anisotropic  ty  any  joint 
set  having  a  preferred  orientation.  The  off-diagonal  terms  of  the  compliance 
matrix  of  an  intact  rock  are  the  same  as  those  of  a  mass  of  the  same  rock 
containing  two  sets  of  orthogonal  joints.  It  is,  therefore,  incorrect  to 
assume  that  Poisson':;  ratio  of  the  jointed  mass  and  the  rock  material  are 
equal.  Table  1  lists  the  particular  expressions  from  which  the  elastic 
modulus  of  the  continuum  anisotropic  rock  mass  may  be  obtained. 
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FIG. A, 2a1  MAXIMUM  PRINCIPAL  STRESSES  (pSl)  IN  A  QUARTER  BLOCK  OF  THE 
ROCK  MASS  SUBJECTED  TO  AN  OVERALL  SHEAR  STRESS  OF  519  PS  I 
(POSITIVE  VALUES  INDICATE  TENSILE  STRESSES), 
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FIG. A. 2b  MAXIMUM  PRINCIPAL  STRESSES  (PSI )  IN  A  QUARTER  BLOCK  OF  THE 
ROCK  MASS  SUBJECTED  TO  AN  OVERALL  HORIZONTAL  STRESS  OF 

-1850  PSI. 
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Table  1 


V  Moduli  of 
>s>r*ock  mass 

Joint  sets 

E,  =  RF,  E 

1  1  r 

E2  =  RF2Er 

V1 

vs 

1 

°12 

Single  joint 
set  normal  to 
the  axis  2 

1 

i 

1  +  Er 

vr 

1+1 

°r  %> 

Orthogonal 
joint  sets 
normal  to 
axes  1 
and  2 

1+  Er 

1*  Er 

VpRPi 

vrRP2 

1+1 

Vm 

S2KN2 

Gr  S2KT2 

Orthogonal 
joint  sets 
with  stag¬ 
gered  cross 
joints  nor¬ 
mal  to  axis 

1 

1  ♦  ¥r 

1  +  Er 

V*! 

i 

1 

1  +  bti  +  1 

S1«N1 

S2KN2 

Gr  S1KT1  G2^J 

Note  -  The  modulus  !i^  and  Poisson’s  ratio  in  the  direction  normal  to  the 
axes  1  and  2  are  equal  to  Er  and  respectively  where,  suffix  1  and 


2  denote  the  axes  of  anisotropy  which  are  noimal  to  the  joint  sets  ] 
and  2; 

IF  =  reduction  factor  defined  as  the  ratio  of  mass  modulus  to  the 
modulus  of  the  intact  rock  material, 

E  =  mass  modulus, 

S  =  joint  spacing, 

E^,  G  ,v  =  elastic  modulus,  shear  modulus  and  Poisson’s  ratio  of  the  Intact 
rock  material, 
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K.^1  KT1  and  K^,  =  normal  and  shear  stiffness  of  joint  set  1  and  2 

respectively , 

and  EL,-^  =  stress  concentration  factors  for  normal  and  shear  stresses 
along  cross  joint  set  1  which  may  be  computed  from  the 
theoretically  derived  relations  for  a  given  joint  offset  S. 


Abstract  of  Publication  15  -  "Continuum  characterization  of  a  rock  mass  - 
Part  II:  Finite  element  analyses"  by  B.  SINGH 

Finite  element  analysis  is  used  to  compare  the  displacements  field  given 
by  the  anisotropic  continuum  model  of  a  jointed  mass  with  that  obtained  from 
tne  discrete  joint  model.  The  models  are  described  in  part  I  of  the  paper. 
Two  basic  problems  were  considered: 

(1)  A  vertical  line  load  acting  on  a  semi-infinite  rock  mass  in  which 
joints  were  oriented  at  0°,  30°,  60°  and  90°  with  the  horizontal; 

(2)  Internal  pressurization  of  a  circular  opening  in  a  rock  mass  con¬ 
taining  orthogonal  joints. 

'hie  analyses  were  made  for  both  continuous  and  staggered  joints.  The 
results  have  been  expressed  In  dimensionless  form  and  may,  therefore,  be  used 
for  other  values  of  joint  stiffness  and  spacing  (see  Figures  A. 3,  A. 4  and  A.b) 
'iho.ie  computations  reveal  excellent  agreement  between  the  finite  element 
predictions  of  the  joint  model  and  the  continuum  model,  except  In  the 
region  of  steep  stress  gradients  near  the  loaded  area.  Tills  conclusion  is 
of  significant  importance  because  a  very  substantial  number  of  nodes  are 
saved  in  doing  away  with  the  joint  elements  in  the  continuum  analysis.  Much 
larger  areas  can  then  be  realistically  simulated  by  the  finite  element  analysi 

Further,  a  simple  method  has  been  proposed  to  predict  the  average  de¬ 
flections  below  a  loaded  area  for  joints  having  arbitrary  orientation  in  all 
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fig.a-3  comparison  of  vertical  deflections  u  computed  bv  finite 

ELEMENT  ANALYSES  USING  THE  JOINT  MODEL  AND  THE  CONTINUUM 
MODEL 
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FIG, A. 5  COMPARISON  OF  RADIAL  DEFLECTIONS  U  COMPUTED  BY  THE  FINITE  ELEMENT 
ANALYSES  USING  THE  JOINT  MODEL  AND  THE  CONTINUUM  MODEL, 


tluve  dimensions.  Tlie  method  requires  the  calculation  of*  the  Increase  In 
.••train  energy'  due  to  the  presence  of  Joints  in  the  rock.  "Three  cases,  of 
various  joint  orientations  and  different  locations  of  the  load  with  respect 
the  joint,  have  been  analysed.  It  was  found  that  the  proposed  method  correctly 
predicts  the  trend  of  the  finite  element  results. 

it  is  felt  that  the  method  is  sufficiently  accurate  to  help  in  idealizing 
.an  actual  tliree  dimensional  rock  mass  in  a  two  dimensional  finite  element 
model.  Thus,  ttie  three  dimensional  deflections  are  first  computed;  the 
equivalent  two  dimensional  finite  element  model  is  then  found  by  trial  and 
error,  being  that model  for  which  the  deflections  correspond  with  those 
Detained  in  the  three  dimensional  case. 

Abstract  of  Publication  16  -  "Reliability  of  dllatometer  and  plate  load 
tests  in  the  determination  of  the  modulus  of  deformation  of  a  jointed  rock 
mss"  by  B.  SINGH 

A  search  of  the  literature  and  enquiries  of  design  engineers  revealed 
that  the  data  available  on  the  results  of  field  plate  jacking  and  borehole 
Jilntometer  tests  is.  very  limited.  It  is  quite  inadequate  for  the  proper 
assessment  of  the  reliability  of  these  test  procedures. 

A  computer  pm-yam  was  therefore  devised  to  simulate  a  jointed  rock 
mass.  Test  data  was  then  'developed'  by  conducting  (i.e.  by  mathematical 
simulation)  a  series  of  plate  load  tests  at  various  surface  location;',  on  the 
ji.ii.nti.-d  mass.  Similarly  dllatometer  tests  were  simulated  and  test  'data' 
obtained  from  the  computer  model  of  the  jointed  mass. 
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FIG.  A. 6  MODULUS  REDUCTION  FACTORS  OBTAINED  BY  FINITE  ELEMENT 
ANALYSES  OF  A  STRIP  LOAD  ACTING  AT  DIFFERENT  PLACES 
ON  THE  SURFACE  OF  ROCK  MASS  CONTAINING  INCLINED  JOINTS. 


FIG.  A. 7  RESULTS  OF  MODULUS  REDUCTION  FACTOR  (INDICATED  BY 

COMPUTER  SIMULATED  TESTS)  MEASURED  BY  A  DILATOMETER 
IN  A  ROCK  MASS  HAVING  RANDOM  JOINT  CHARACTERISTICS. 
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The  deflection  of  the  instrument  was  computed  for  random  spacing, 
orientation,  and  stiffness  of  the  joints,  and  random  location  of  the 
instrument,  These  results  were  then  analyzed  to  determine  the  minimun 
number  cf  tests  required,  in  terms  of  the  average  joint  frequency  and  the 
modulus  reduction  .actor  [i.e.  the  ratio  of  in-situ  modulus  to  intact  specimen 
m  lulus.  I,  in  order  that  the  error  in  predicting  the  mass  modulus,  would  be 
Jess,  th.in  30:7 ,  with  a  confidence  interval  of  90#  (see  Figure  6  on  page  32  in 
main  text). 

Lt  was  concluded  that  the  optimum  orientation  of  the  dilataneter  to 
tne  joints  is  in  the  range  of  10°  to  30°  (see  Figure  A.7).  The  number 
of  tests  required  increases  rapidly  with  decreasing  modulus  reduction  factor. 

Similar  computations  were  made  for  plate  load  tests.  (See  Figure  A. 6). 

It  was  , noted  that  the  mass  modulus  is  not  affected  significantly  by  the 
sine  of  the  plate. 

Statistical  moments  of  the  test  results  were  found  to  be  of  considerable 
n,'l;  in  interpreting  the  test  data.  Tiius,  the  standard  deviation  of  the 
le  flection  was  found  to  be  almost  independent  of  the  size  of  the  loaded  area, 
ili is  enables  a  distinction  to  be  made  between  scatter  due  to  erroneous 
measurement  and  that  resulting  from  the  stochastic  nature  of  the  rock  mass. 
Further,  the  third  moment  of  the  deflections  increases  sharply  as  the 
f  the  loaded  area  decreases, .  'lids,  may  enable  the  i vaults  of  small 
scale  *>'sts  to  ie  correlated  with  tiius.o  oJ’  large  scale  tests. 

ip.p-w  Interpretation  of  the  scatter  in  small  scale  test  data  c;ui 
cens  i  terui  l.v  enlianc-'  the  value  of  the  results. 


